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Abstract 

This paper is devoted to the study of the emptiness formation probability 
r (to) of the open XXZ chain. We derive a closed form for r (to) at A = 1/2 when 
the boundary field vanishes. Moreover we obtain its leading asymptotics for an 
arbitrary boundary field at the free fermion point. Finally, we compute the first 
term of the asymptotics of In (r (to)) in the whole massless regime — 1 < A < 1. 



1 Laboratoire de Physique, ENS Lyon et CNRS, France, karol.kozlowski@ens-lyon.fr 



1 



1 Introduction 



For already a few decades, spin chains trigger a lot of interest as many relevant 
physical quantities can be computed. For instance the spectrum ( [1], [2]) of the spin- 
1/2 XXZ chain as well as its correlation functions ([3], [4]) are known. The open XXZ 
spin-1/2 chain with diagonal boundary conditions: 

M-l 

n =Ys a >™+l + < a m+l + A ~ !) + <A h - + °M h+, (1.1) 

m=l 

is another example of an integrable spin chain Its spectrum was first obtained by 
Alcaraz et al. [5] in the framework of coordinate Bethe ansatz. One year later, 
Sklyanin [6] developed an appropriate scheme to apply the algebraic Bethe ansatz 
to models with boundaries. In the mid nineties, the Kyoto group obtained an inte- 
gral representation for the elementary blocks of a half-infinite massive (A > 1) XXZ 
spin-1/2 chain with a diagonal boundary field [7]. Recently, the Lyon group derived 
formulas for the elementary blocks of the finite open chain XXZ with diagonal bound- 
ary fields and for an arbitrary anisotropy parameter [8]. 

The simplest possible correlation function of a spin-1/2 chain is the so-called empti- 
ness formation probability (EFP) r (m) first introduced in [9]. This quantity can be 
understood as the probability of observing a ferromagnetic string of length m starting 
from the first site of the chain. The EFP can be computed to the end in the case of 
a periodic XXZ spin-1/2 chain at A = 1/2. Indeed this exact expression has been 
conjectured in [10] and proven in [11]. The free fermion point of the bulk model is also 
quite particular as the EFP can be represented as a Toeplitz determinant of a smooth 
function on an arc. Hence its asymptotics can be studied [12] by using the strong 
Szego limit theorem for circular arcs [13]. One can also derive the leading asymptotics 
of lnr (m) of the bulk massless model using a saddle point method [14]. 

The aim of this paper is to obtain corresponding results for the open XXZ spin 
1/2 chain. The symmetrized integral representation for the EFP derived in [15] will 
be the starting point of our considerations. We first show that r (m) can be computed 
to the end in the case of an open spin-1/2 chains subject to a zero boundary field 
at A = 1/2. After setting some notations in Section [21 we will present this result 
in Section [3] We then derive the leading asymptotics of the EFP at the free fermion 
point in Section [U This is done by exploiting the almost Hankel determinant structure 
of the EFP at A = 0. We analyse the asymptotic behavior of such determinants in 
Appendix EJ This analysis takes advantage of the asymptotics of Toeplitz matrices 
with Fischer-Hartwig type singularities [16] as well as of the uniform asymptotics 
of orthogonal polynomials of modified Jacobi type [17]. In the last section we use 
the saddle point approximation borrowed from the bulk [14] to derive the leading 
asymptotics of lnr (m) in the massless regime. 
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2 The Emptiness formation probability. 

The authors of [15] showed that, in the case of a half-infinite chain, the emptiness 
formation probability admits a multiple integral representation whose integrand is 
a symmetric function of the integration variables. This way of representing r (m) 
enables the separation of variables in the m-fold integral defining r (m) at A = 1/2. 
It also allows to express r (m) as an almost Hankel determinant at A = 0. Finally, it 
is suited for the derivation the leading asymptotics of In [r (m)] in the massless regime. 

We recall the standard parameterizations of the massless regime 

A = cos(C) ,CG]0;vr[ and h- = - sin C cot £_ , |_ <E ] -tt/2 ; vr/2] 

as well as the symmetrized integral representation for the emptiness formation prob- 
ability [15]: 



r(m) = {E&K..d m) 



— [ 

n ml J 



22 / 

d m A det m [£ (Aj,Cfc)] det m [l/h(\j,£ k )} 



n* 2 

C(h) i<j 



Q— 7 ~ ^ (2-1) 



5 (X jk , i() a {X jk , i() k=i a (A fc , i|- + </2) 

The contours of this multiple integral depend on the boundary magnetic field h- : 

C (h) = R U r± |t* (€- + C/2) } whenever - C/2 < |_ < and C (/i) = R otherwise. 

Here T± stands for a small circle skimmed through in the anticlockwise/clockwise 
direction. Furthermore, we agree that 

5 (x, y) = sinh (x + y) sinh (x — y) c (x, y) = cosh (x + y) cosh (x — y) 

Xij = Aj — Aj A^ = Aj + Aj 



and introduce the functions: 

p (A - - P (A + sinh (vrA/C) sinh «/C) 



sinh (2£) . 

C sinh (2^) c 



(2.2) 
(2.3) 



h(X,0 = S (A + e,</2)s(A-e,</2) (2.4) 

Lastly, we remind that p (A) = — — — — solves the Lieb equation: 

C cosh (7TA/C) 

sin ( /* . sin 2£ 



vrs(A,iC/2) 



{X) + J dpK{X-p)p{p) ; K(A) = 2 ^ A 2C <} (2.5) 



3 



3 C = 7r/3, an exact result. 

When C = 7r/3, the integrand of eq. (|2.ip can be further simplified due to the 
duplication formula: 



sinh (3x) = 4sinh (x) sinh {x + in/ 3) sinh (x — in/3) 
Indeed, in this special point one recasts r (m) as 



(3.1) 



T (?7l) 



(3/4) 



m(m+l) 



i<j 



d m A 



to! {An) 



C(h) 



n 

k=l 



sinh (3A fc ) sinh (2A fc ) s £ k , i£ 



s ( Xk,ii- + ivr/6 



xdet r 



1 




1 


det m 


.c(Aj,Cfe). 





(3.2) 



A partial homogeneous limit has already been performed in the latter formula. 

The symmetry of the integrand enables us to replace the Cauchy determinant 
det m [1/c (Aj, by to! times the product of its diagonal elements. This yields the 
sought separation of variables. In the case of a zero boundary field ( ie. = in/2), 
the evaluation of the integrals leads to 



r (to) 



m(m+l) 



(-2)' 



det m [g (£j,Zk)] 



fl sinh (2£ fc ) sinh (60 ,0,* 2 &'6 



fc=l 



where 



g(x,y) 



sinh (x + y) /2 



j<k 

sinh (x — y) /2 



(3.3) 



(3.4) 



sinh3(x + y)/2 sinh 3 (x — y) /2 
It happens that it is possible to evaluate the homogeneous limit of (13. 3|) . Indeed 
we have the following equality of limits: 

det m [g(Cj,Ck)] __ r det m [g(xj,y k )} 



lim ■ 



= lim 

ue k uu"-e k ) 2 5 u ■<!,. w u 



k=l 



HI 



(3.5) 



k=l i<j 

We can send the x'a and the y's to zero in an arbitrary way when computing the 
limit in the RHS of (|3.5p . In particular, we can choose xi = ai and yj = a (m + j), 
and sent a to 0. Such a homogeneous limit is the a — ► limit of det m [Z7], where 
J7 G 7W m (C) is given by 

sinh a (i + j + m) /2 sinh a (j — i + to) /2 



sinh ft {i + j + to) /2 sinh (3 (j — i + m) /2 
Actually Kuperberg [18] evaluated det m [U] for all a and j3: 



(3 = 3a. 



(3.6) 



det m [U] 



2m 2m+l 

J! 2 sinh /J (j - i) /2 II 2 sinh (a + (3(J - t))/2 

2|j 

m 

Q 4 sinh /3(m + j - i) /2 sinh (3(m + j + i) /2 

i,3 



(3.7) 
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Hence: 



lim m detm }^ )] 2 = lim^^M (3.8) 



n«n(g-g) n(* 2 -i 2 ) n (* 2 -j 2 ) 

fc=i «<j j<j <j<j 

1< m+l 

This allows to obtain a closed formula for the emptiness formation probability in the 
homogeneous limit: 

2m 2m+l 

nu-i) n i+3(j-i) 

,(„.)-- S££ - M (3 9) 

" ' (-2) m (2m)! m m 2m ' 

' X[(m + j - i) (m + j + i) ]\ (j 2 - i 2 ) J] (j 2 - i 2 ) 

i,j i<j <i<j 

m+l 

Or in terms of factorials: 

-w=^n(3,-i)n^Ti5! (^n^rw ■ (310) 

It is now a matter of standard asymptotic analysis to extract the large m behavior 
of (CTH : 

T (ro)=(^Y m^(l + o(l)) , (3.11) 



where the constant C is expressed in terms of the Euler Gamma function T(z), the 
Euler constant 7 and of the Riemann Zeta function^ (z). 

c - iW exp G c ' ( - I) + 3l (1+ ^) (312) 

+ OO 

r ^ <r* /2 - cosh t/3 -cosh t/6 + 2sinhf/3 + sinht/6 5 



1 - e~* V t 2 6t 72 







We stress that the qualitative behavior of r (m) differs from the bulk one by the 
presence of an exponential factor e cm in the asymptotics . This seems to be a general 
feature of the boundary model as a similar behavior appears at the free fermion point. 
Moreover the gaussian decay is twice faster than in the bulk. 

4 The leading asymptotics of r (m) at A = 

The integral representation for the EFP also admits a separation of variables at 
A = (ie C = 7r /2). More precisely specifying eq. tj2.lh to ( = tt/2 we get 
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2m(m+l) 2 s j n 2m ( £ 
t M = jc—yii -det m [T] 

The entries of the matrix T are defined by an integral: 

sinh 2 (2A) 



cosh fc+ ' (2A) cosh (2A) + sin ( 2£ 



dA 



(4.1) 



(4.2) 



We perform the change variables z = 2/ cosh (2 A) — 1 in the above integral. Then 
t (to) reads 



t (to) 



and 



det m [R] 
(-27r/i_)" 



where Rjk = e* 2 lim 



J+k-2 



0+1 z — a 



cj(z)dz , (4.3) 



]-l-ie;l-ie[ljr + (o) h- G ] 1 ; +oo [ 
] — 1 — ie ; 1 — ie [ /i_ G ] — oo ; 1 [ 



(4.4) 



(z) = y/(z-l) (3 + z) and a = /i_ + /C 1 - 1 



X + X 1 



In order to apply the results of Appendix[A]we recall the Wiener-Hopf factorization 

(4.5) 



of v3~+ cos 9 

y/3 + cos I 

Specializing eq. ()A.35|) to the determinant representation for the EFP, we get 

(h-XT 



Tim) ~ C — 2 — < 
v ' 2 m V to / * 



'1 + x" 1 ^ 1 



1 + X 



-i 



. m / — 1 , —1 

X \ X +u 



h- < 1 



h_ > 1 



(4.6) 



h-l V 1 + x- 1 
The constant C is expressed in terms of w and of the Barnes G-function 

.., _ ,'1 + U~ l Y 7TlG(l/2) 



1-u" 1 / Vl + 



(4.7) 



One can see explicitly from the large to asymptotics of r (to) the asymmetry in- 
duced by the boundary field. This asymmetry appears for | h- |> 1. One has 
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T h- ("») / (» + X + ')(l + X-') , , m , , .... 

" V ("xi 1 + !)(! + *;■) (_X+X_) ( ' 

where 

X± ^ X±1 = ± + /C 1 ) - 1 . (4.9) 

Finally, one infers from the asymptotics that r (m) — ► for large m, as it 

should be. Indeed, in such a limit, the first spin is necessarily oriented upwards. 
Conversely, when h- — > +oo, r (m) is not vanishing, also as expected. 

5 The rn — > +00 limit of In (r (m)) /m 2 

Using the saddle point approximation we obtain the leading asymptotics of In (r {rn)). 
We find that the leading terms is gaussian and independent of the value of the bound- 
ary field h- : 

In (r(m)) = Km 2 + o (m 2 ) (5.1) 

The constant K is twice bigger than the one appearing in the leading asymptotics of 
r (m) of the periodic chain [14]. 

In order to implement the asymptotic analysis of In r (m) we first recast r (m) as 



si nh(g-) \ 2m ( it \ m /t^ 2 "' 
(2() 2m \sm(()J U, 



r(m) = ^^n (__) («) Jm (5 .2) 



I m is an m-fold integral 



T 



J d™Ae m MMr) detm [M jfe ] (5.3) 



2 



Vcosh 2 7rA fc /Cy ' 
and the matrix M is obtained by inverting the equation 
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We now reexpress Z m in a form more suited for an asymptotic analysis: 
d m Ae™MMr) detm[M ] +S ^"WC 



0<Ai<-<A ft 



sinh(2A m ) 



im-l\ 



0<Ai<-<A m— 1 
-iA ro =C/2+|_ 



xe 



m 2 ^({A>r- 1 ) 



det r 



M 



5 (A m ,iC/2)5 2 (ttC^AjQ/C) 
(Aj m , iC) s (Ajm, «C) cos]l4 (vrA m /C) 



Here 



1 £-€ 
otherwise 



(5.6) 



(5.7) 



Note that we have used the symmetry in the A's as well as the A — ► —A invariance 
of the integrand in order to replace the integration over R m by 2 m m\ times the inte- 
gration over the ordered domain {0 < Ai < • • • < A m }. We also used the determinant 
structure and the symmetry of the integrand to evaluate the contribution of poles, if 
it exists. 



The leading contributions to the above integrals will be equal to the integrand 
evaluated at the solutions of the saddle point equations. We assume that, in the 
m — > +oo limit, these solutions densify on M + with a density a (A). Then sums can 
be replaced by integrals according to the prescription 



i=i i 



dAa(A) / (A) 



(5i 



It is worth noticing that in this limit, Mj p can be approximated by 

K (A 7P ) - K (X ip ) , 
M jp = S jp + K 3V) V 3V) . 5.9 
ma (Xj) 

This can be seen by replacing the sums by integrals in (15. 5p and using the integral 
equation (|2.5p . Then, Hadamard's inequality 



det m [a,jk] \< I max (| a jk 



m 2 



(5.10) 



ensures that 



lndet m [M,- fc ] 
nm k — - — = nm 

m— >+oo m m— >+oo 



In det r 



M 



j,k 



. 



As a consequence, the functions det m [M] and det T 
leading asymptotics of In (r (m)). 



M 



(5.11) 



cannot contribute to the 
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The density of saddle point roots a (A) satisfies a singular integral equation. 
This equation is obtained by taking the large m limit in the saddle point equations 

1 =0. This singular integral equation reads: 



2 ( ^ tanh ^ - coth (A + i(/2) - coth (A - i(/2)\ = (5.12) 

/27T \ d 
dfia (fi) ~r tarih 7T — coth (A — e/i + iC) coth (A — e/z — i£) 



It is natural to extend the density into an even function on R. This recasts the 
integral equation into a form very close to the integral equation appearing in the bulk 
model [14]: 

tanh - coth (A + i(/2) - coth (A - i(/2) = (5.13) 

V.P. j d/icr (/i) | — tanh tt — — coth (A — fi + ^C) c °th (A — yu — *C)T • 

R 

Eq. (|5.14l) is solved by the Fourier transform. We find 

2coshfcC/2 , m V2co S hvrA/2C 

" (fc)= coshfcC h6nCe " (A)= CcoshvrA/C (5 - 14) 

It is then immediate to evaluate the m — ► +oo limit of ({A}™) at the solutions 
of the saddle point equations: 



]]m ? (I \W =?{<*>) = f dk cosh2 ( fc C/2) sinh [k (tt - C) /2] 
m ^ +oc ^ ^ J-i ; j k cosh (jt0 sinh (K/2) sinh ((fcvr)/2) ' 1 J 

R-iO 

The two multiple integrals appearing in (|5.6|) have thus the same gaussian decay. 
Thus, 

lim i^Ml = 21nf^+^ M (5-16) 

In particular we recover the results at A = and A = 1/2. As already announced 
we obtain twice the corresponding bulk constant. Unfortunately the saddle point 
method cannot catch the exponential behavior of the leading asymptotics of r (m) 
and in particular the dependence on the boundary field. 
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6 Conclusion 



We have obtained the leading asymptotic behavior of the EFP of the open XXZ 
spin-1/2 chain for some particular values of the couplings. The A = 1/2 point seems 
to have, at least at h- = 0, similar symmetry properties to the bulk model. Thus 
it should be possible to evaluate to the end the integrals appearing the generating 
function of a z correlators at A = 1/2 in the open chain. As suggested by the saddle 
point analysis, the presence of the boundary increases the long range gaussian decay 
of r (m) by roughly a square factor with respect to the bulk. It would be interesting 
to derive some formula for the exponential term e cm appearing in the asymptotics 
of In (r(m)) in the massless regime. All the more that this terms should exhibit a 
dependency on the boundary field /i_ . 
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A Appendix 

In this Appendix we derive the leading asymptotics of a modified Hankel determi- 
nant detmR. The entries of R reads 

Rij = e i7ra lim / dz z [ +k ~ 2uJ ^ . (A.l) 

^ n i z - a P> 
P =i 

The integration is carried over a segment and loops encircling some poles of the inte- 
grand: 

% = [-l-ie;l-ie] [j T + (a p ) with £c[1;j]. (A.2) 

In what follows, we assume a p G C \ ] — 1 ; 1 ] and t E whenever a p £ ] — oo ; — 1 ]. 
Moreover we restrict to weight functions of modified Jacobi type: 

u>{z) = {z-lT{l + zfg{z) . 

We also assume g (z) holomorphic in an open neighborhood of 

s 

[-1;1] U and such that g([-l;l])c K + . 

p=i 
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A.l Reduction to a simpler problem 

We first remove the eventual contour integrals from most of the lines of R. For 
j < m — s, a linear combination of Lj, . . . , Lj +S allows to perform the replacement 



z J+k ^ z 3+k JJ (1 _ z J ap) 
p=l 

in the integrand. For j € [ m — s; m — 2] one can only replace 



(A.3) 



] J (i _ z / ap ) 

p=j+2—m+s 



The above transformations lead to 



s 

det m [R] = det m [#] (-a p ] 



s+l—m—p 



p=l 



where 



i? 



/ dzz- ? ' +fc - 2 U( 



3 



1 < j < m — s 



z j+k-2 ^ t z \ 

e io7r li m + f m — s < j < m 

u j j- m + s J — 

* n (^-op) 
P =i 



(A.4) 



(A.5) 



(A.6) 



and uJ(x) = e J7ra u; (a; - z0 + ) 

The second step separates the "pole" part of the determinant from its pure Hankel 
part. Let IL^ (z) be the monoic orthogonal polynomials on [— 1;1] with respect to 
the weight U(x). The reconstruction of these polynomials in the first m — s lines as 
well as in all the columns makes the first m — s lines diagonal. This procedure splits 



det r 



R 



into a product of two determinants. Namely, 



det T 



R 



= det m _<j [H] det s 



(A.7) 



Where 



K 



j z k+j - 2 uJ(z) dz 
-l 

^ Ilfc+m— s— 1 



jk 



lim 



/ 



uj (z) dz 



n (z- ap) 
p=i 



(A.8) 
(A.9) 
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A. 2 The Hankel Determinant 



The asymptotics of the Hankel determinant can be obtained via a procedure es- 
tablished by Basor and Ehrhardt. They built a mapping between Hankel and Toeplitz 
determinants [19]. This mapping allowed them to infer the large size behavior of the 
Hankel determinant from the known asymptotics of the Toeplitz determinant. We 
recall their result in 

Proposition A.l Let T and H be respectively a Toeplitz and a Hankel matrix 

TV 

— [ c(9) e-^-Ofl 'd9 (A.IO) 
2vr J 



- kl 



— TV 

1 

k+l-2 



H M = ZJ(x)x k+i - z dx (A.ll) 



-i 



such that uj G L 1 [ — 1 ; 1 ] and 

c{9) = isgn(9)uo (cos 9) - vr < 9 < it , (A. 12) 

then 



det m [H] = 2- m ^ m - l \ m ^/det 2m [T] (A.13) 

The function c : 9 i— > isgn (9) uJ (cos 9) is of degenerate Fischer-Hartwig type [16] 
as it admits two maximal decompositions into the canonical Fischer-Hartwig product: 

c (9) = 2- a -?h (cos 9) w^ 1/2 (e w ) u> ,«,-i/2 (e*) 
-2- a ^h (cos 9) u> WlA _i/ 2 (e ie ) u> , a) i/ 2 (e") 
The function W0 rjTj)J (e* 61 ) reads 

^0 r ,r, CT (e ie ) = (2 - 2 cos (0 - # r )) r e **V-*-0r) 9 r < 9 < 9 r + 2ir (A.14) 

The assumptions made on g (z) guarantee the existence of a Wiener-Hopf decom- 
position: 

g(cos9) = C\g}g + (e i6 )g + (e- ie ) , (A.15) 



where C [g] and g + (t) satisfy the constraints: 

C[g]eM g + (t)=exv\^t n [g} \ ; [<?] n e R (A.16) 



, n=l 



x If a or /3 = —1/2 then c(6 l ) is not degenerate. In this case one should drop the ( — l) m factor in 
(| A.17[) . Moreover the result doesn't rely on a conjecture at these points as asymptotics of Toeplitz 
determinants with such Fischer-Hartwig symbols are known [16]. 
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The leading order asymptotics of H thus follow from the generalized Fischer- 
Hartwig conjecture raised by Basor and Tracy [20]: 



detm lTu (e)1 = (my g m w ^-'c ± lit + „ (A17) 

where 

TT G(l + t+l/2)G(l + t-l/2) , _.. X>[<^ , 

/<*,/? = II G(l + 2t) ^[5]=e 1 . (A.18) 

Here G is the Barnes G function. 

The asymptotic behavior of the Hankel determinant reads: 



det m _ s [H] = 2- {m ~^ 



( w) \ 0, n am V gw (1 + 0(1)) 

v - ' m4^(-l)^(l) 



(A.19) 



A.3 Asymptotics of det, [i^ m )] 

Computing the poles and taking the e — > + limit yields 



1 z i+m-«-i nfc+ro _ a _ 1 ( z )^ 

J 

]1 (Z - Op) 

p=l 



^ = /- A k+m - s - l[z) u(z)dz (A.20) 



+ W Q7r fc+m ' s ' l(af) u; (a t ) 1 B (t) 



n a* 

p=l 

¥* 

Here lg stands for the characteristic function of E. The asymptotic analysis of (IA.21|) 
is based on the uniform asymptotic estimates for the monoic orthogonal polynomials 
with respect to the modified Jacobi weight U [17]. Let K be compact in C \ [— 1 ; 1] 
and K compact in ] ; n [. Then : 

n n (cosfl) = / J P °° 2l/2 n = cos(fn+ " + f + 1 )e + *(e)\e n (A.21) 



yju (cos 9) sin (9) 
= « = ^i.lxl>X. (A.22, 
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One has e n , e' n =^1 + y — ) j uniformly in € K and respectively : e K. Moreover 
have introduced the functions: 



i 



,W = 2^ -2U + a| (A ' 23) 



= VHs+^^-^Ta+X- 1 )" (A.24) 

D ~ = J^^vlS (A - 25) 

Recall that we have decomposed z into z = (x + X -1 ) /2, | X |> 1 and Wiener- 
Hopf factorized g. It follows immediately from (|A.22|) that the leading order of the 
"pole part" in (|A.21|) reads 

m— s+j— Itt z' n /ir l m—s+j—1 i /r,\m—s+k—l/2 

a t n jfc+m _ s _i (a t ) C [5] a t (Xt/ 2 ) 



3 2 a +' 3 //-, -2\ 



a; (a t ) 



p=l 



/ x (Xt ~ !)" (xt + l)f / (t . nR s 

x 5+ (Xt) " e n ( A -26) 

n a t P 

P =i 

¥* 



We have used the identity 

Ut-lT(xt + lf (a t -l) a {a t + if 



(A.27) 



2a+ " (i-XtT^ + Xt 1 )" 

We now derive the asymptotics of the integral term in (|A.21 j) . This amounts to 
the study of the asymptotics n,q — > +00 , n — q + j > 1, of the sequence 

1 

^,» = / f ""^ 57(z)dz (A.28) 

-1 U(z- a p ) 
p=i 

Applying the formula for the asymptotics of the orthogonal polynomials (1A.21|) we 

get: 



_ Doo /" „ COS 9 6* (COS 9) , „ N ' 3 4- 1 , ^"-i- 1 

T = — 2L. / dfl (1 + cos 6»)2 + 4 1- cos 0)2+3 







^ l{n+ l±^ )e+me) + ^(^A^)^^ ^ (A 2Q) 
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We recast the integral as one over ] — it ; it [ and then, using the analycity of the 
integrand, we shift the contour to 9 € ] — it + i0 + ; it + i0 + [. Since 

(C ° S ^" 1} "\ (l-cos0)V^ ]-tt + *0+ ;*)+[ (A - 30) 

we can absorb the factors e T ^( a+ ^) . Moreover 2"> (cos ± l) 7 = (l + ±e if? ) 7 G 
] -it + i0 + ; vr + i0 + [. Finally, we obtain 



c lg] jT(i+^(i-^r i ~"^( e .)^ (A . 31) 

-7r+io+ n ( c ° s # - °p) 
P =i 

The above integral can be interpreted as a contour integral over T (0, 1 — + ). 

j c\g\ f dC (i + C) m/2 (i - C) a+1/2 n (nrnf C-' + C Y r 

x 9,n - 2 n+a+/3+l f iC 3 y+VSH I ^ I n 

r(0,i-0+) I! (C + C- 1 -2a p )/2 
p=i 

The integrand has no pole at C, = 0; indeed it behaves as £-? -1+n- P, £ — > 0, and we 
have imposed n + j — p— 1 > 0. The only poles inside of the contour are in the points 

C = xt 1 ,t e 1 1 ; j 1 where 2a t = xt + xi 1 and I xt |> i- So 

^^t^ + ^) (l -^ (1+ f'^ ^ (A, 2 ) 

4=1 x? (xr 1 - x t ) n a* 

P =i 



T - _ _ 



Thus has the asymptotic behavior 



^ 2 *=i 

iT^' -1 (i + o(i)) 



(x^xr^a-x^^i + xr 1 )"} 



\/i - x t 2 n <% 



p=i 



One can simplify the asymptotics of det s [i<C( m )] by performing the linear combi- 
nation of lines: 
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u 



3-1 j-l+t 

L '-Y.—j — 
- 1 n a* 

P=t+i 



L t 



j = 2...s 



(A.33) 



Then 



det. 



s(s + l) 
2 2 



1 



in 



m—s+t—l 



(l + o(l)) 



(A.34) 



xdet s 



*e i7r V ( Xt ) xT' s+k - 1 (Xt - IT (1 + Xtf 1e (t) 

-9+ (Xt 1 ) Xt 



Note that the terms coming from the integration over ] — 1 ; 1 [ are exponentially sub- 
leading with respect to the ones coming from the residues. Nevertheless we keep them 
in the final formula as we could have t E. For instance, in the limit E = we 
recover the expected formula (1A.19|) for the asymptotics of Hankel determinants with 

_ 3 

a weight ul/ Y\ ( z — a p)- 
P =i 

Summing up all the results, we get the formula for the leading asymptotics of 
det m [R] 



det m [R] 



5+1 



irC[g] 



in 



-2 — 



4 (-1)^(1) n 2a fcp 

p>k 



fafiE [g] 



det. 



\j n(i-x t - 2 ) 

9 + {XtlXl {l-XtTi l + Xt l Y (A.35) 
ie i ™g + ( X t) XT'^ 1 (Xt ~ If (1 + Xtf 1e (*)1 (1 + o (1)) 
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